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ABSTRACT 


The  effect  of  various  jamming  waveforms  on  conventional 
binary  incoherent  digital  receivers  was  analyzed  in  terms  of 
resulting  receiver  performance  (i.e.,  receiver  probability 
of  error ) . 

Probability  density  functions  associated  with  the  test 
statistic  generated  by  incoherent  receivers  under  the  influ¬ 
ence  of  noise  and  jamming  have  been  obtained. 

Due  to  the  complexity  of  the  mathematical  expressions 
specifying  receiver  probability  of  error  in  closed  form,  no 
attempt  has  been  made  to  obtain  absolute  optimum  jamming 
waveforms  operating  against  binary  incoherent  receivers. 
Therefore  near  optimum  jammer  signals  were  proposed, 
studied,  and  evaluated. 

The  effect  of  a  varying  threshold  on  receiver  perform¬ 
ance  was  investigated  and  a  jamming  strategy  involving  use 
of  an  FM  jammer  was  considered,  and  its  effect  evaluated. 

Graphical  results  are  presented  that  highlight  the  math¬ 
ematical  results  obtained. 
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I  .  INTRODUCTION 

A.  COHERENT  CORRELATION  RECEIVER 

The  optimum  receiver  which  will  detect  the  presence  of  a 
signal  or  discriminate  between  two  different  signals  in  the 
presence  of  additive  white  Gaussian  noise,  is  the  well-known 
coherent  correlator  receiver  structure  which  has  two  equiva¬ 
lent  forms,  as  shown  in  Figure  6.1  [Ref.  2]. 

The  binary  communication  problem  is  modeled  using 
hypothesis  testing  theoretic  principles  in  which  one  of  two 
signals,  S| (t)  or  S0(t),  is  received  in  the  time  interval 
[0,T].  Due  to  the  presence  of  the  noise,  the  observable 
signal  r(t)  takes  on  one  of  the  two  following  forms 

H,  :  r(t)  =  S,  (t)  +  n(t)  CKt  4T 


or 


.  r(  t )  -  S9(  t )  +  n(  t )  0  1 1  ^T 


where  the  noise  n(t)  is  assumed  to  be  a  sample  function  of  a 
stationary  white  Gaussian  process  having  zero  mean  and  power 
spectral  density  (PSD)  level  N0/2  watts/Hz. 

The  optimum  receiver  (in  the  sense  of  minimum  prob¬ 
ability  error,  Pe  )  generates  the  test  statistic  G,  where  G 
is  given  by 


(r,Cj ) 


I 

2- 


II SJI  -  Is.  II 


or  equivalently 
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and  compares  it  to  a  threshold  in  order  to  decide  whether 
signal  Si(t)  or  S0(t)  has  been  transmitted  in  the  time 
interval  [0,T],  The  inner  product  notation  (•,•)  implies 

(x,y)  =  f  x(t)y(t)  dt 


O 
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We  can  interpret  E  as  the  average  energy  per  bit  and  Q  as 
the  normalised  signal  cross  correlation.  Using  the 

following  definition  of  the  error  function 

f  i 

erfc^Cx)  -  J- 

Equation  1.1  can  be  written  as 


From  the  above  equation  it  can  be  noted  that  as  (1-P)E/N„ 
increases,  the  error  probability  decreases.  For  fixed  E/N0 
therefore,  the  optimum  signals  choice  is  one  for  which  P  = 
-1.  This  occurs  whenever  S,(t)  =  -S0(t).  A  system  employing 
this  choice  of  signal  is  known  as  an  optimum  antipodal 
signaling  binary  communication  system. 


B.  JAMMING  OF  COHERENT  RECEIVERS 
1.  General 

The  objective  of  this  section  is  to  summarize  previously 
derived  results  on  the  performance  of  the  coherent  binary 
(.optimum)  receiver  of  Figure  6.1  operating  in  a  jamming 
"ent .  That  is,  the  transmitted  signal  is  interfered 
by  the  presence  of  a  jamming  waveform  as  well  as  additive 
white  Gaussian  noise  [Ref.  1],  The  signal  appearing  at  the 
front-end  of  the  receiver  can  be  mathematically  modeled  by 

r(t)  =  S,(t)  +  n(t)  +  n.  ( t )  04t6T  , 

J  "  (1.2) 

r(t)  :  E0(t)  +  n(t)  +  n.(t)  0£t6;T 

J 

where  n.(t)  is  a  jammer  waveform  modeled  as  deterministic, 

0 

yet  unknown  to  the  receiver. 
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The  resulting  average  probability  of  error  is  given  by 


=  pR'P(M°i  +  P~  ?(H|)  771 

=  KM.)  /  )  <•*  »  pen, )  / cg| h  ,  )  eg 

f)  -co 


where  71  is  the  threshold  with  which  G  is  compared  in  order 
to  decide  which  is  the  true  hypothesis.  Observe  that 

0)  t@(0)  d£> 

The  first  integral  in  the  above  expression  for  Pe  is 
(probability  of  false  alarm).  It  can  be  expressed  as 
f o 1  lows 


where 

)  Is(  aV)  d V 

is  a  well-known  tabulated  function  called  the  Marcum 
Q- function.  The  second  integral  in  the  above  expression  for 
Pe  is  PH  (probability  of  miss).  It  can  be  expressed  as 
f o 1  lows 


(2.3) 


25 
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where  S9  is  used  in  place  of  sin(Wc  t  +  0  ).  It  can  also 
be  seen  that 


Var  r  E  |(n,S)  | 


M.T  ,  sin  2WeT 


-<|hJ 


and  similarly 


Var  <  Y 


sin  2W. 


For  convenience  we  assume  WcT  =  n7C,  where  n  is  an  integer. 
Thus  the  sine  terms  above  vanish  resulting  in  var(X  I H#  *  0  > 
=  var(Y  |  hj,  .  d  )  =  N0  T/4  =  (J  2  .  If  this  assumption  is  not 
made,  an  additional  term  results.  However,  if  7L  /Wc  <<  T, 
the  additional  sine  term  is  small  and  may  be  neglected. 
Furthermore  the  covariance 


E  |[x  -  E  {x|H.,0}]  •  [  Y  -  E  {y|H;,0J  ]  I  H.,0 
=  E  j(n,S)  (n,C)|  ‘  "  1 


N0  T  1  -  cos  2W/r 


i  =  0,1 


provided  the  assumptions  on  Wc  hold.  This  implies  that  for 
any  given  value  of  ©  ,  both  X  and  Y  are  uncorrelated 
Gaussian  random  variables  and  therefore  statistically  inde¬ 
pendent.  The  density  function  of  G2  conditioned  on  the 
phase  ©  is  noncentral  Chi-Squared  distributed,  and  is  given 
by 


23 


r(t )  =  n( c)  +  n.  ( t ) 

j 


0  £T 


where  n  .  ( t )  is  the  jammer  waveform  present  during  the  time 
int e rva 1  [ 0  ,  T] . 

In  the  absence  of  n.(t),  the  optimum  receiver  for  the 
binary  ASK  problem  is  well-known.  Its  derivation  is  well 
documented  in  the  statistical  detection  theory  literature 
[Ref.  2].  The  receiver  structure  is  shown  in  Figure  6.2 

In  this  section,  the  effect  of  n.(t)  on  this  receiver  is 

o 

analyzed  by  evaluating  the  resulting  Pe,  under  the  assump¬ 
tion  that  n.(t)  is  a  deterministic  jammer  waveform,  however 
unknown  to  the  receiver  itself.  Receiver  performance 
requires  determination  of  the  statistics  of  either  Gl  or  G, 
where  GJ  is  the  output  of  the  quadrature  detector  and  is 
given  by 


where 


X  =  £r(t)  sin  W-t  dt  *  (r,.l) 


Y  =  i  (t) 


(r,0 


Provided  that  the  random  variable  £p  is  fixed  to  some  value 
$,  X  and  Y  are  conditional  Gaussian  random  variables  with 


E  f X|H,  ,ej  =  (AS8 ,S)  ♦  (a.  ,S) 


rii  =  ,C)  +  (n.  , 


‘J’c) 
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II.  ANALYSIS  OF  JAMMING  ON  ASK 

A.  GENERAL 

In  this  chapter,  the  effect  of  a  deterministic  jammer 
waveform  on  the  performance  of  an  incoherent  ASK  receiver 
will  be  investigated. 

The  signal  at  the  front-end  of  the  receiver  is  given  by 
Equation  1.8  with  the  modification  that  under  either 
hypothesis,  a  jamming  waveform  n;(t)  is  present  during  the 
time  interval  [0,T], 

For  the  coherent  ASK  receiver,  the  optimum  energy 
constrained  jammer  waveform  could  be  obtained  by  using 
Equation  1.5  .  However,  the  complexity  of  the  mathematical 
expression  for  probability  of  error  in  the  noncoherent  case 
makes  it  very  difficult  if  not  impossible  to  derive  the 
optimum  energy  constrained  jammer  waveicrm  in  closed  form  as 
will  be  seen  in  the  sequel. 

A  reasonable  postulation  is  that  the  optimum  jamming 
waveform  for  the  coherent  ASK  receiver  can  act  as  a  good 
jammer  for  the  incoherent  ASK  receiver  also.  Thus,  such 
near  optimum  jammer  signals  are  studied  and  evaluated  in 
terms  of  their  effect  on  the  performance  of  the  noncoherent 
ASK  receiver. 

B.  ANALYSIS  WITH  NEAR  OPTIMUM  JAMMER 

Analysis  of  the  incoherent  receiver  starts  from  the 
mathematical  model  of  the  receiver  front-end  input  waveform 
r ( t  )  given  by  either 

!!,  :  rft )  =  A  sin  (Wct  +  if)  )  +n(t)  +n.(t)  O^t^T 
or 
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where  we  assume  that  the  signal  amplitudes  even  though 
random,  remain  constant  over  a  T  -  second  interval.  The 
noise  is  modeled  as  in  the  previous  analysis,  and  the  signal 
amplitudes  have  (a  priori)  probability  density  functions 
given  by 

fA(a)  =  -71  exp  (-  ^  ~  0 

*  h.  o  -  r-o 

and 

b 

f„(b)  =  —77.  exp  (-  -0  n)  b  >  0 

B  1  o 

where  =  E(AJ)  =  E{B2}. 

The  results  associated  with  ASK  incoherent 
receivers  can  be  applied  toward  obtaining  the  optimum 
incoherent  FSK  receiver.  This  is  worked  out  in  Reference  2. 
The  resulting  receiver  structure  is  shown  in  Figure  6.5  . 
Individual  channel  matched  filters  shown  in  Figure  6.5  can 
be  essentially  replaced  with  tuned  bandpass  filters  having 
center  frequencies  corresponding  to  the  "mark"  and  "space" 
frequencies  W,  and  W<>. 


an  envelope  detector  and  a  sampler  as  shown  in  Figure  6.4 
[Ref.  2]. 


3.  FSK 

For  binary  frequency  shift  keying  (FSK),  the 
received  signals  at  the  front  end  of  the  receiver  in  the 
presence  of  noise  only  can  be  mathematically  modeled  by 
either 

H,  :  r(t)  =  A  sin  (W,t  +0))  +  n(t)  QAt^T 

or 

H0  :  r(t)  =  A  sin  (W0t  +  (j)  )  +  ntt)  0^-t^T 

where  the  additive  noise  is  again  assumed  to  be  zero  mean, 
white  Gaussian,  with  power  spectral  density  level  N0/2 
watts/Hz.  Usually  the  frequencies  W(  and  WQ  differ  substan¬ 
tially.  The  phases  (Q)  and  <£)  are  statistically  independent 
random  variables,  uniformly  distributed  over  the  interval 
[0 , 2 7t  ]  • 

Insofar  as  the  signal  amplitudes  is  concerned,  two 
cases  are  distinguished.  In  the  first  case,  the  information 
bearing  signal  amplitudes  are  known  and  equal.  In  the 
second  case,  the  information  bearing  signal  amplitudes  are 
treated  as  independent  random  variables  that  are  Rayleigh 
distributed  due  to  multipath  propagation.  For  this  second 
case,  the  following  hypothesis  testing  theoretic  model  is 
used 

H  |  :  r(t)  =  A  sin  (W,  t  +  0  )  +  n(t)  0  4  t  AT 


and 


H 0  :  r(t)  =  B  sin  (W0t  +  (f)  )  +  n(t)  QAt^T 
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2.  ASK  (On  2  Off  keying) 

In  the  case  of  incoherent  ASK  (On  -  Off  keying),  the 
signals  at  the  front-end  of  the  receiver  in  the  presence  of 
noise  only  can  be  mathematicaly  modeled  by  either 

(1  8 ) 

H,  :  r(t)  =  A  sin  (W£t  +  0)  +  n(t)  O^t^T 
or 

Hc  :  r(t)  =  n(t)  0  irt 


The  amplitude  A,  the  frequency  Wc  and  the  time  of  arrival 
are  assumed  to  be  known  except  that  the  phase  <0  is  modeled 
as  a  random  variable  having  an  a  priori  density  function 


0  £  0  1  Z7TL 

otherwise 


The  additive  noise  is  again  assumed  to  be  zero  mean,  white 
Gaussian,  with  power  spectral  density  level  N0/2  watts/Hz. 

In  the  absence  of  a  jamming  waveform  n.(t),  the 
optimum  receiver  for  decoding  the  binary  information  (which 
is  transmitted  via  the  use  of  either  a  sinusoid  that  is 
incoherently  received,  or  no  signal  at  all)  in  each  interval 
of  duration  T  sec,  is  the  well-known  quadrature  receiver. 
Its  structure  is  shown  in  Figure  6.2  [Ref.  2]. 

There  is  an  alternate  form  of  the  quadrature 
receiver  obtained  by  replacing  the  correlators  with  matched 
filters  having  an  impulse  response  given  by  h(t)  =  sin  Wc(T 
-  t)  and  h(t)  =  cos  Wc(T  -  t),  0  <  t  <  T,  and  sampling  the 
outputs  at  t  =  T  (Figure  6.3)  [Ref.  2]. 

Still  another  important  alternate  form  of  the  quad¬ 
rature  receiver  is  an  incoherent  matched  filter  followed  by 
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and 


O':  1  -p  = 


(A,  -  Aa* 

X.  X- 

A,  +  Ac 


The  'break  point'  for  ASK  occurs  at  JSR  =  0^/2.  Since  qC  < 
1,  in  terms  of  'break  point'  efficiency,  PSK  is  highest  with 
a  'break  point'  occuring  at  a  JSR  of  1.  FSK  is  next  highest 
in  efficiency  with  a  'break  point'  occuring  at  a  JSR  of  1/2, 
and  ASK  is  lowest  in  efficiency  with  a  'break  point’ 
occuring  at  JSR  <  1/2. 


C.  INCOHERENT  RECEIVERS 
1 .  General 

In  the  coherent  systems  considered,  the  information 
bearing  signals  were  assumed  to  be  known  exactly  at  the 
receiver . 

In  noncoherent  systems,  the  phase  of  the  carrier 
signals  is  not  available  at  the  receiver  so  that  the  phase 
is  treated  as  a  random  variable  uniformly  distributed  over 
[0 , 2  7L] •  As  such,  we  may  expect  the  performance  of  an 
incoherent  system  to  be  degraded  in  comparison  to  the 
performance  of  the  corresponding  coherent  system.  However, 
because  of  their  simplicity,  incoherent  systems  are  widely 
used  in  many  applications. 
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with  the  constraint  that  (W|  -  W0)T  =  1  71  ;  (W,  ♦  W0)T  = 

m7L ;  1  and  m  integers.  Finally  for  ASK  modulation 


S,(t)  =  A,S(t)  ;  S0(t)  =  AeS(t)  CUt^T 


where  we  assume  that  ||S||  <  co  and  for  convenience 
>  A0 .  The  optimum  jammer  waveforms  and  resulting 
performances  can  be  obtained  using  Equations  1.5 
Thus  the  optimum  jammer  for  PSK  is  given  by 


n  .(t) 

j 


0  it  6  T 


Since  P  =  - 1  for  PSK,  we  have 
CD 


P  =  — 

t?  2 


1  e-xV2 


V27C 


J2SUR  (1  -y.JSR) 


dx 


For  FSK  modulation,  the  optimum  jammer  is  given 


n.Ct)  = 
j 


7=r  sin-^CW,  -  Wa  )t  cos  (W,  +  W0  )t 

J  JT 


Using  the  assumptions  on  W(  and  W0 ,  the  value  of  P 
Thus  Equation  1.7  becomes 


1 


?  = 


xX/2 


-y-sr/R  ci  +  y?jsR  ) 


IjTjz  Q 


-  xx/2 


dx 


P7L  ^  dx  + 

:',r'R  fl  -  /2JSR  ) 

Finally,  for  ASK  modulation,  the  optimum  jammer  is 

S(t) 

n.Ct)  = 

j 


and  the  probability  of  error  becomes 


,  that  A( 
receiver 
and  1.7. 


is  zero. 


given  by 


Analysis  of  Equation  1.7  indicates  the  fact  that  for  a  value 
of  JSR  beyond  a  so-called  'break  point',  that  is  JSR  > 
(1-/0  )/2,  an  increasing  SNR  causes  Pe  to  become  worse,  i.e., 
it  increases. 

3.  Optimum  Jammer  Waveforms  for  P SK ,  FSK  and  ASK 

The  effectiveness  of  a  deterministic  jammer  waveform  on 
various  modulation  techniques  will  now  be  presented  for  PSK 
modulation,  where 

S,(t)  =  A  sinWct  S0(t)  =  A  sin(Wct  +  7C)  O^t^T 

with  the  constraint  that  WcT  =  n 7c  ;  n  an  integer.  Also, 
results  for  FSK  modulation  will  be  shown,  where 


3,(t)  =  A  sin  W,t  ;  S0(t)  =  A  sin  W0t 


0  At 
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|d|  =  |(n;,Sd)  <  flnj.flsj 


with  equality  if  n.(t)  is  proportional  to  Sj ( t ) .  Defining 

j~  "  'J 

|  |  n .  |  |  =  /  (2  , where  P*.  is  the  jammer  energy,  the  term 

j  *  1  j  j 

|  d  |  -  -  >  oo  implies  that  Pn-  -  ->  oo  when  ||Sd||  <  oo  .  Since 

■J 

it  is  not  possible  to  have  infinite  jammer  energy,  one  must 
constrain  n.(t)  such  that  P*.  is  finite.  From  the 
Cauchy -Schwarz  inequality,  Equation  1.4  can  be  made  into  an 
equality  if  - 


n j  C  t )  =  k.SjCt) 


where  K  is  a  constant  of  proportionality. 

K  must  be  set  to  the  value  J  P„7  /  |  |  Sj|  |  . 

J 

ized  by  setting 


If  |  |n.  |  |  =  J?Z  t 

V  J 

Thus  | d |  is  maxim- 


"j(t>  -J^T  sd(t)/  ||Sj|| 


(1.5) 


and  from  the  above  discussion,  this  results  in  Pe  being 
maximized  also.  Substituting  Equation  1.5  into  the  prob¬ 
ability  of  error  expression  (Equation  1.3)  yields 


It  is  possible  to  put  Equation  1.6  in  a  more  meaningful  form 
by  defining 

E  ^ 

r  J-  "  •  signal  to  noise  ratio 


where 


%  S 


;  d  =  (n.  ,S.) 


From  the  jammer  standpoint,  the  optimum  jammer  waveform  must 
be  chosen  in  such  a  way  that  it  maximizes  the  receiver  prob¬ 
ability  of  error.  A  determination  of  the  optimum  jammer 
waveform  can  be  made  by  first  evaluating  the  derivative  of 
Pe  with  respect  to  d,  which  is  the  cross  correlation  between 
the  jammer  waveform  and  the  signal  difference  Based 
on  Equation  1.3,  we  can  show  that 


*  d‘]/2 


Due  to  the  behavior  of  the  hyperbolic  sine  function,  we 
observe  that 


cKi 


>0  ,  d  >  0 

=  o  ,  d  =  0 

<0  ,  d  <0 


421 r 


[  -  Kj  M  *  "//  2  ]  > 


for  all  values  of  d.  Since  Pe  is  a  monotonic  function  of  d, 
it  is  apparent  that  making  d  as  large  as  possible  in  magni¬ 
tude  results  in  the  largest  possible  increase  in  Pe.  In  the 
limit  as  |  d  |  -  -  >  cp  ,  it  is  easily  seen  that  Pe  -->  1/2. 
However,  from  the  Cauchy  -  Schwarz  inequality 


\ 


2.  proper  Optimization 

The  effect  of  a  jammer  waveform  on  the  optimum  cohe¬ 
rent  receiver  is  now  presented  in  terms  of  its  impact  on  the 
receiver  probability  of  error. 

In  order  to  determine  the  performance  of  the 
receiver,  we  analyze  the  test  statistic  G,  which  is  given  by 


G  =  <r,Sj)  *4r|Kf  -  ||S.f] 


with  r(t)  now  given  by  Equation  1.2  .  Due  to  the  assumption 
of  a  Gaussian  zero  mean  noise  having  PSD  level  N0/2 
watts/Hz,  it  is  simple  to  show  that  G  is  a  conditionally 
Gaussian  random  variable  with  mean  value 

E  (g|SxJ  =  (S;,Sj)  +  (n.  ,Sj)  *  VLIN|  *  IIs.  II  ]  ’  1  =  °’1 


and  variance 


■Ms.} 


IT  lpJ 


i  =  0,1 


With  the  mean  and  variance  known,  the  resulting  conditional 
density  function  of  G  is  given  by 


VglSi)  = 


exp  - 


'27C*Ne||^ 


(S;,Sj)  -  (r..,Sj)  - 

"  2-%-  s 


From  the  general  probability  of  error  expression  (Equation 
1.1)  and  an  approriate  change  of  variables  in  the  integrals, 
when  P ( H ,  )  =  P(H0)  =  1/2  the  expression  for  Pe  becomes 


as 


m 


fcO 

frf 

f  (0) 

.  (gli 

i  ,8)  dg 

0 

IJ 

G  1  9 

loo 

0 

- 

d0 


where  in  the  second  equality  the  order  of  the  integration 
has  been  changed.  The  inner  integral  of  Equation  2.3  can  be 
expressed  in  term  of  the  Marcum  Q-function  as 


^  1  exp  (-  )  I0(iSGl 


1  a 1 


)  dg 


m. 


V  exp  (-  — ---  —  )  L(  aV)  dV 


=  1  -  Q  (  a  ,  n  /  a  ) 

,  r~ 

where  =  \J$/(y  <  Therefore  the  probability  of  error  can  be 
written  as 


=  P(H0)  Q  (  —  ,  n  /a  ) 
e  o 


(2.4) 


+  P  ( H  ,  ) 


CzK 


1  - 


ys- 


-'0 


o  (  ,  n  /a  )  d  e 

^  o 


where  the  dependence  on  0  is  imbedded  in  the  term  / S  and 
the  threshold  TJ  for  the  jammer  absent  case  is  obtained  as 
the  solution  to  the  equation 


-  A  T  /  2  M<j 

g  I0(  2An/W„)  =  R 


R  = 


P(H0) 
P  (  H  ,  ) 


which  can  be  equivalently  put  in  the  form 

W 


e 


A2'T/2Ma 


I, 


T/N,  (  H  /  o  )  |  =  r 


'1 


The  I0  ( •  )  function  used  here  and  also  previously  used  in 
conjunction  with  the  development  of  Pe  is  the  modified 


26 


Bessel  function  of  the  first  kind.  The  appropriate  setting 
of  the  threshold  is  an  important  issue  for  both  coherent  as 
well  as  noncoherent  receivers.  One  could  use  the  threshold 
setting  that  would  be  derived  from  the  analysis  of  receivers 
operating  in  additive  white  Gaussian  noise  interference 
only.  This  approach  can  be  quite  unsatisfactory  as  demon¬ 
strated  in  Reference  1.  A  better  approach  would  be  to 
obtain  Pe  as  a  function  of  the  threshold,  and  then  search 
for  the  threshold  that  minimizes  Pe.  While  this  approach  is 
intuitively  appealing,  it  is  often  mathematically  intrac¬ 
table.  While  the  threshold  issue  is  not  addressed  in  this 
particular  section,  it  will  be  discussed  in  more  detail  in 
the  sequel  and  simulation  results  will  be  presented. 

If  the  definition  of  average  signal  energy  previously 
introduced  is  used,  we  have  E  =  A*T/4  which  is  reduced  in 
half  in  comparison  to  the  binary  signal  transmission  case, 
due  to  the  fact  that  the  information  bearing  signal(s)  do 
not  have  equal  energy.  In  order  to  afford  comparisons  with 
the  coherent  receiver  case,  we  will  implicitely  boost  the 
value  of  signal  amplitude  A,  to  obtain  E  =  AzT/2  in  order  to 
have  agreement  with  previous  cases  insofar  as  signal  ener¬ 
gies  is  concerned.  Thus  the  threshold  determination  equation 
now  becomes 


P(H„)  (2.5) 

r(H, ) 


2 2 TIP  (  n  /a  )  1 


1 


If  we  are 
maximize  Pe 
3P 
3S  = 


to  find  the  optimum  jammer  waveform 
,  an  attempt  must  be  made  to  solve 


0  and 


3P 
3S/  = 


0 


so  as 


to 


Unfortunately  the  resultant  equations  are  mathematically 
involved  and  do  not  appear  readily  solvable  for  n.(t).  It 
seems  however  that  a  'good'  jammer  waveform  can  be  postu¬ 
lated  based  on  the  results  obtained  for  coherent  ASK.  It  was 
found  for  that  case  that  the  optimum  n.(t)  is  a  tone  at  the 

v) 

carrier  frequency.  Thus  the  following  jammer  waveform  can  be 
used  as  a  potential  near  optimum  jammer,  namely 


0  i  t 


(2.6) 


Observe  that  with  this  choice,  |  |  n  •  |  |  2  =  P-n.  .  The  prob- 

ability  of  error  Pe  can  now  be  determined  using  the 

threshold  setting  equation  (Eqn.  2.5)  and  the  previously 
derived  expressions  for  S  and  S'.  The  effect  of  the  near 
optimum  jammer  waveform  on  the  receiver  (i.e.,  incoherent 
receiver  performance)  can  be  analyzed  by  evaluating  Pe  as  a 
function  of  n.(t)  using  Equation  2.6  .  Note  that  in 

Equation  2.6  the  jammer  energy  is  E,,.  .  It  can  be  shown  that 

J 

(nJ  »s)  =  [Kff  sin1  Wct 

(nj,C)  =JT sin  Wct  cos  Wct  dt  =  0 


Then  from  Equations  2.1  and  2.2,  the  probability  of  error  Pe 
given  by  Equation  2.4  can  be  expressed  in  terms  of  SNR 
( s ignal- to-noise  ratio)  and  JSR  ( j ammer- t o - s ignal  ratio) 
using  the  fact  that 


28 


7lT  '  2  SNR  (1  +  2  J  'lSK  cosQ  +  J5R) 

=  2  SNR -JSR 

and  the  fact  that  ^  &  can  be  obtained  from  Equation  2.5. 

For  equally  likely  hypotheses,  the  probability  of 
receiver  error  (Equation  2.4)  can  be  expressed  by 


pe  =  +  Q  (  y2SMR  J3R  .  7^/(T ) 

j>sm  (1  +  2v'T5R  CCS0+  JSR  ) 


(2.7) 


From  Equation  2.7,  Pe  greater  than  0.5  can  occur  if 

Q  (a  ,8  Q  (a'  ,  8  >  d0 


(2.8) 


where 


CX  =  J2SUR  (1  +  2/TSR  ccsO  +  JSR) 

Ct  -  J 2SNR~JSR  and  -  Tj  /q- 

Since  the  periodic  sinusoidal  function  imbedded  in  c(  '  is 
integrated  for  one  period,  the  condition  oC  -  '  might 

cause  the  condition  of  Equation  2.8  to  be  satisfied  for  a 
fixed  value  of  ^  .  Thus  there  is  a  possibility  of  obtaining 
Pe  greater  than  0.5  for  a  JSR  value  beyond  the  critical 
value  of  0.25  which  can  be  obtained  from  the  condition 
ol  -  oC  '  ,  that  is  , 

/Tgnr-jsr  =  (1  +27ynT"co:.e  *■  jgk) 

/ICR  COS0  =  -  J_ 


(2.9) 


Receiver  performance  in  the  presence  of  a  jammer  is 
expressed  in  terms  of  signal-to-noise  ratio  (SNR)  and 
j ammer- to- s iganl  ratio  (JSR).  A  rectangular  pulse  of  dura¬ 
tion  T  seconds  has  amplitude  spectrum  AT  sine  Tf,  and  B  = 
2/T  is  a  rough  measure  of  its  bandwidth.  Thus  the  expression 
of  the  form 

A  T/2  A  /2  A  /z  (2  10) 

No  N.  2  PSD  -  B 

2.  'T 

can  be  interpreted  as  the  ratio  of  signal  power  to  noise 
power  (SNR)  in  the  signal  bandwidth  [Ref.  5].  The  term  JSR 
can  be  expressed  by 

P«.  _  jammer  energy  (2.11) 

AiT/2  signal  energy 


where  P^.  represents  the  jammer  energy  defined  before. 

J 

C.  VARIABLE  THRESHOLDING  EFFECT 

It  is  apparent  that  when  a  fixed  threshold  value  is  used 
by  the  receiver,  the  effect  of  the  jammer  waveform  having 
energy  such  that  the  resulting  JSR  value  is  above  a  certain 
level,  is  such  that  the  receiver  may  be  rendered  inoperable. 

It  could  therefore  be  suggested  that  setting  the 
threshold  based  on  the  value  P(H|)  or  P(H0)  only  (see  Eqn. 
2.5),  may  not  be  desireable.  Values  of  threshold  other  than 
some  fixed  value  may  result  in  improved  receiver  perform¬ 
ance,  in  other  words,  reduce  the  jamming  effect. 

Recall  that  the  threshold  setting  equation  is  of  the 

form 

-SNR  _ 

2  I0(  /''—P  a)  =  R 
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where  o(  =  ^  /(J  •  we  now  can  attempt  to  reduce  Pe  by  an 
appropriate  choice  of  threshold  ^  .  By  varying  the  value  of 
R  (instead  of  using  R  =  1)  we  can  obtain  this  threshold 
which  denote  q(  '  from  the  expression 

SNR 

i;'(R6  ) 

where  I0'l(X)  is  the  inverse  modified  Bessel  function  of  the 
first  kind.  Using  an  approximation  to  I0 (X)  which  is  given 
by 


for  x  1 


Equation  2.12  can  be  rewritten  as 


+  j2w(2SNR)  +  4 SNR  +  jtf\  R 


(2.13) 


From  Equation  2.13  ,  it  can  be  recognized  that  when  SNR  is 
large,  variation  of  the  value  of  R  does  not  significantly 
affect  the  threshold  value  tyi  '  because  taking  the  logarithm 
of  R  reduces  the  effect  of  R  further.  Thus  the  large  value 
of  SNR  suppresses  the  effect  of  variation  of  the  term  J2n  R 
on  .  This  limited  variable  thresholding  effect  on  the 
receiver  which  is  under  significant  jamming  is  analyzed  in 
graphical  form  in  Chapter  5  for  various  values  of  R. 
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III.  ANALYSIS  OF  JAMMING  ON  FSK 


A.  GENERAL 

For  binary  incoherent  FSK  with  a  jammer  present,  the 
received  signals  under  the  two  hypotheses  are  either 

H,  :  r(t)  =  A  sin  (W,t  +  0  )  +  n(t)  +  n.(t)  0  it 

j  ~ 

or 

H 0  :  r(t)  =  A  sin  (VI 0t  +  (£>)  +  n(t)  +  n.(t)  0  It  i T 

By  separating  the  frequencies  W,  and  W0  sufficiently,  we  can 
form  signals  that  are  orthogonal,  have  equal  energy,  and 
have  the  same  advantage  of  ease  of  generation. 

The  modified  FSK  receiver  structure  which  is  capable  of 
varying  the  output  of  the  each  envelope  detector  (which  is 
followed  by  a  multiplier)  is  diagramed  in  Figure  6.6  .  The 
optimum  receiver  for  the  case  where  no  jammer  is  present  can 
be  derived  from  statistical  decision  theory  and  is  shown  in 
Figure  6.5  .  This  receiver  can  be  obtained  by  setting  o<(  = 

1/2  in  the  modified  FSK  receiver  shown  in  Figure  6.6  . 

In  practice,  incoherent  FSK  is  widely  used  because  of 
its  simple  receiver  structure,  its  small  performance  penalty 
due  to  lack  of  phase  coherence  and  its  more  efficient  use  of 
signal  energies  in  comparison  to  incoherent  ASK.  In  addi¬ 
tion,  we  are  not  faced  with  the  difficulty  imposed  by  a 
threshold  that  must  change  with  SNR  as  is  the  case  with 
incoherent  ASK.  Thus,  a  receiver  which  is  known  to  be 
optimum  for  incoherent  FSK  transmission,  has  been  modified 
by  including  some  channel  weighting.  This  has  been  done  in 
order  to  be  able  to  determine  whether  or  not  such  channel 
weighting  can  reduce  the  effect  of  the  jammer. 


32 


This  chapter  is  devoted  to  investigating  the  performance 
of  the  modified  incoherent  FSK  receiver  in  the  presence  of 
jamming  and  additive  white  Gaussian  noise.  By  letting  ~ 
1/2  we  can  as  a  byproduct  obtain  the  performance  of  the 
conventional  incoherent  FSK  receiver  of  Figure  6.5  in  the 
presence  of  jamming  and  additive  white  Gaussian  noise. 


B.  ANALYSIS  WITH  NEAR  OPTIMUM  JAMMER 

The  modified  receiver  performance  can  be  obtained  by 
introducing  a  null  hypothesis  (no  signal)  as  a  dummy 
hypothesis  and  by  following  the  same  reasoning  as  in  the 
analysis  of  incoherent  ASK  presented  in  the  previous 
chapter. 

The  receiver  function  is  to  compare  the  envelopes  at  the 
output  of  each  channel  once  every  T  seconds  and  decide  in 
favor  of  the  larger  of  the  two  envelopes  (Figure  6.6).  For 
the  purpose  of  analysis,  let  us  first  assume  that  a  'mark' 
signal  has  been  transmitted,  that  is,  the  hypothesis  H,  is 
assumed  to  be  true.  An  error  is  committed  if  Vc  exceeds  Vf  . 
An  error  is  also  committed  if  Vt  is  larger  than  V0  when  a 
'space'  signal  has  been  transmitted,  that  is,  the  hypothesis 
H0  is  assumed  to  be  true  [Ref.  3]. 

Let  Pe,  denote  the  probability  of  the  first  type  of 
error  described  above,  which  is  expressed  as  Prob.(V0  >  V(  | 
Hj).  Under  the  assumption  that  a  'mark'  signal  has  been 
sent,  the  output  qf  of  one  of  the  envelope  detectors  is 
given  by 


qf  = 


where 


X,  =  rft)  sin  W,t  dt  ;  Y,  =  /  r(t)  cos  W.t  dt 
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Observe  that  X(  and  Y|  conditioned  on  the  phase  and  either 
of  the  two  hypotheses  are  Gaussian  random  variables  with 


|H,,0}=  / A  sin  (W  ,t  +  0  )  sin  VI ,  t  dt 

+  fTn.(t)  sinW.tdt  =  (S,  ,S),  +  (n.,S)t 

>0  *  J 


and 


^  {  A  I  Wi»©|=  [a  sin  (W,t  +  Q  )  cos  W,t  dt 

fr  £> 

+  I  nj(t)  cos  W,t  dt  =  CS  ,C),  +  (n.,C), 


where  Sj  represents  the  function  A  sin  (W, t  *  6  ),  S  repre¬ 
sents  the  function  sin  W,  t  and  C  represents  the  function  cos 
W,t.  Likewise,  assuming  again  that  W,T  =  n;L,  n  an  integer, 
and  that  n(t)  is  zero  mean  white  Gaussian  noise  with  PSD 
level  N0/2  watts/Hz,  we  obtain 


Var  jx,  |  H,  ,©}  =  Var  {yJh,,©}  = 
Furthermore,  it  can  be  shown  that 


h0T 


|[.X,  -  E  |x,  (  H,  ,  0}  j  [Y|  -  E  |  Y,  (  H,  ,6}] 


H  i  ,0 


=  0 


so  that  the  conditional  r.v.'s  X,  and  Y|  are  uncorrelated 
and  therefore  independent.  The  sum  involving  random  vari¬ 
ables  X,  and  Y , ,  and  producing  q*,  will  result  in  a  non¬ 
central  Chi-Squared  distribution  so  that 


1 


f  0t|9  (q  H|  }  =  2oa'  f" 


q,  +  s.,  j 
TTT)  lo(^) 


q, 


whe  re 


s„  =  rz{y,|H(  ,ej  +  Ea  {  y,  I  h , , e) 


[(S,  ,S),  ♦  cnj,s),]  4-  [fs,  ,o,  4  cnj,o; 


34 


where  the  dependence  on  0  is  imbedded  in  the  term  Sn  .  We 
now  need  to  obtain  the  statistics  of  the  output  of  the 
multiplier  following  the  upper  envelope  detector.  (See 
Figure  6.6).  That  is,  the  probability  density  function 
fV|  (V,  |  H()  needs  to  be  derived.  From  standard  transforma¬ 
tion  theory,  using  the  relation  V,  =  2(1  -  &L  )q( ,  it  can  be 
shown  that 


On  the  other  hand,  the  output  of  the  lower  envelope 
detector  when  H  |  is  assumed  to  be  the  true  hypothesis,  is 
given  by 
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,,  z-  2 
/.  +  Y 

o  L  0 


where 


•"•<>  =  J  sin  Wat  dt  ;  Y 0  -  J  r(t)  cos  W„t  dt 

Following  a  similiar  procedure  as  used  above,  it  can  be 
shown  that 


e|y0|H,  ,ej  =  j  A  sin  (W,  t  +  0)  cos  W0t  dt 


+  n  .  ( t )  cosDW  t  dt  =  (S,,C)  +  ( n .  ,  C ) 


E{xo|H(,eJ  =  I  A  sin  (W,  t  +  0)  sin  WDt  dt 

+  f  nJ(t)  s in 0W  t  dt  =  (S,,S)0  +  (n  .  ,S)fl 

Here,  S  and  C  have  the  same  meaning  previously  defined 
except  that  the  subscript  "0”  outside  the  inner  products 
implies  that  we  should  interpret  S  as  sin  W0t  and  C  as  cos 
Wct.  It  can  also  be  demonstrated  that 

Van  |  Xo|H,,0}=  Var  {y0|H,,qJ  = 

and  also  that 

E  {  [  x0  -  E  |  y.J  H(,qJ  ]  [  Y„  -  E  [  Yo|H,,eJ]|H,,0j=  0 

so  that  the  conditional  r.v.’s  X0  and  Y0  are  uncorrelated, 
hence  independent.  Thus  similiar  to  Equation  3.1,  the 
expression  for  the  conditional  density  function  of  qc 
becomes 
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where 


q0^° 


E1|xo|H(,0j  +  iJ 

[<3,  ,S)0  +  (nj,S) 


H,  ,©j 

[(S.  A 


+  (n.  ,C)0J 


z 


and  dependence  on  0  now  is  imbedded  in  the  term  S  oj . 
Applying  again  the  random  variable  transformation  V0  = 
2  (/q0 ,  the  conditional  density  function  of  V 0  is  given  by 


f  Vq(  V0|  H  ,  ) 


1  p*  Vo 

Urn  / — exp  (- 

Jo  a  K 


(3.3) 


Va  _ 

I0(  2 


In  order  to  compute  the  probability  of  error,  we  can  now 
use  the  previous  expressions  for  the  conditional  probability 
density  functions  which  are  derived  assuming  a  'mark'  signal 
has  been  sent.  That  is,  for  a  given  value  of  V,,  an  error 
is  made  if  V0  >  V|.  Thus  the  average  error  probability  is 
found  by  averaging  the  conditional  error  probability  given 
by 


-30 


PC/,)  =  f  (V0|H,,V.)  dV„ 


over  all  Vi .  That  is, 


(3.4) 


i 


i 


i 
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(Vo  H,  )  dV0  fV|  (V.  I  H  .  )  dv, 


Substituting  for  the  density  functions  in  Equation  3.4,  from 
Equations  3.2  and  3.3,  we  obtain 


(271.  z 

I  l  J  IFexp  (  -  fa2" )  *  S"  V 


L  -'o 


(3.5) 


Ip  (  22  d’/pjde]  fV|(V,|H, 


)  dV, 


where  the  order  of  the  integration  has  been  changed  for 
computational  convenience.  In  the  above  equation  the  inner¬ 
most  integral  can  be  expressed  in  terms  of  the  Marcum 
Q- function  as  follows 


Vo  /  Vp  )  -  2o_  /P— 

~T  exp  (-  ?  i  '  So<  )  I0  (  2a  Vbo*)  dv0 

-V.  0  2ax  o*- 


=  q  (  -2-  ) 

x  a  /  a  a 

Then,  Equation  3.5  for  Pe  becomes 


<  H  fv.<v'|w  av' 


fix  nr— 

1  o  (  ,  — —  X  )  d0 

—  '  a  *  a 

*7  '  n 


(3.6) 


\±n 


!  XP  (  - 


_  4. 

c  *  (~r~)  )  I0C  x  — )  d 0 


where  in  the  second  equality  the  change  of  variable  x  = 

[  ,  has  been  made . 

2  (  1  -  a  )  a 

From  the  orthogonality  property  of  the  signal  pair  used 
(which  is  obtained  by  assuming  sufficient  separation  between 
two  frequencies  and  that  W  ,  as  well  as  W0  are  large),  we 
have 

J  sin  (W ,t  +  0 )  sin  W at  dt  =  0 

so  that  the  term  S  o)  is  independent  of  @  .  Therefore 
Equation  3.6  can  be  rewritten  in  the  following  form 


CD  - 

Q  (  x)  x  exP  (“ 


({El)]  (3-7> 


To  (x 


d? 


where  the  order  of  integration  has  been  changed. 
Furthermore ,  using  the  following  formula  involving  an  inte¬ 
gral  of  a  Marcum  Q-function  (Ref.  A] 


r>  R 

—  )  on  exp  (- 


p.; 


V  A)  dR, 


the  inner  integral  in  Equation  3.7  can  be  simplified  in  such 
a  way  that  Fe,  becomes 


(3.3) 


la(y.  — )  dxl  de  "^e>:p(-  rTl)  dA 


,  p*  r  rco  r~  x 

Q(^,x)  x  exp (  — 


n exp  |  - 


Ao/s — ^  -2 
— rVin  )  -ri. 


Ia(x§  dA  |  dx  d8 


where  the  last  equality  has  been  obtained  by  interchanging 
the  order  of  integration.  This  equation  does  not  appear  to 
be  readily  s imp  1 i f iab le .  If  we  take  a  closer  look  at  the 

innermost  integral  in  Equation  4.1,  that  integral  might  be 
evaluated  using 
rv 

/  u  - 1  ,  o  o 

t  T  (ct  t)  exp(-p  t  )  dt 


—  )(y~) 


2pap(v  +  l) 


exp( 


where  (•)  is  the  modified  Bessel  function  of  the  first 
kind  of  order  v,  p'  ( • )  is  the  Gamma  function,  and 
,F,  (X,,X2;Xj)  is  the  Confluent  Hypergeometric  function.  In 
the  special  case  of  X,  =0,  the  Confluent  Hypergoemetric 
function  ,F,  ( 0 , Xz ; )  =  1  so  that  the  integral  in  question 
can  be  simplified  in  this  case.  However,  since  the  term  Sit 
includes  the  integration  factor  A,  we  may  be  able  to  calcu¬ 
late  and  express  Equation  4.1  in  a  simpler  form  using  the 
formula  above  for  the  case  in  which  the  variable  A  could  be 
separated  out  in  the  term  ^  Su  for  a  given  jammer  waveform. 

The  same  arguments  apply  to  obtaining  the  probability  of 
error  under  the  hypotheses  H0 .  Thus  Pe  can  be  shown  to  be 
given  by 


Jo 

bd!£,x)x  exp  (  -  —  ) 


(4.2) 


.2.  exp 


_ z 

)  -  A 


I0  (x 


)  dA  )  dx  jf 
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H,  :  r(t)  =  A  sin  (W,  t  +0)  +  n(t)  +  n.(t)  0 

J 

or 

H0  :  rCt)  =  B  sin  (Wat  +  ^  +  n(t)  +  n.(t)  O^t^T 

where  the  amplitudes  A  and  B  are  independent  identically 
distributed  random  variables  having  density  functions 


2 


fA(a)  - 

“f*  e  x  p  ( - 

aZO 

rx> 

*  *0 

fB(b)  = 

"^2  exp  ( - 

-o 

-£l  s 

J 

*-  •  ‘  0 

b>0 

with  E{ A}  =  E(B)  =  f  JS.  A0  and  E(A*>  =  E{B*)  =  2A* 

V 


Note  that  FSK  with  fading  is  the  same  as  what  was  previ¬ 
ously  analyzed  (FSK  without  fading)  except  that  now  the 
signal  amplitudes  are  random  varibles.  Therefore  for-  a 
fixed  value  of  A,  the  error  probability  when  a  "mark”  signal 
has  been  transmitted  is  the  same  as  that  of  the  noncoherent 
(nonfading)  FSK  case.  That  is,  we  can  use  the  results  of 
the  preceding  chapter  and  in  particular  make  use  of  Equation 
3.7,  to  obtain  (for  the  special  case  of  =  1/2) 


P  (A) 

e, 


x 


10  ( x  7  )  d  x  d  0 


where  the  dependence  on  A  is  due  to  the  fact  that  the  term 
Sii  defined  previously  depends  on  A.  Then,  Pe,  becomes 

lOO  * 


pe(A)  PCA)  JA 


(4.1) 


•  r-  *  (m 

/  v :  ci  \  ,  x  +  \  a  K 

A  — —  ,x)x  exp(-  — - ) 
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IV.  ANALYSIS  OF  FSK  IN  THE  PRESENCE  OF  JAMMING  AND  FADING 

A.  GENERAL 

In  certain  propagation  media,  the  received  signal  trans¬ 
mitted  via  a  free  space  channel  is  often  subject  to  fading, 
a  phenomenon  caused  by  multipath  propagation  and  the  equiva¬ 
lent  addition  of  random  phasors.  The  vector  sum  signal  will 
have  an  envelope  which  changes  with  time  resulting  in  an 
effect  known  as  fading. 

The  fading  signal  model  often  utilized  is  that  of  the 
nonselective ,  slow  fading,  Rayleigh  distributed  signal 
amplitude  where  it  is  assumed  that  the  amplitude,  while 
random,  remains  constant  over  the  time  interval  [0,T].  Thus 
the  received  envelope  is  now  random  with  probability  density 
function 

9  2- 

fA(a)  =  ~~  exp ( -  )  a  1>  0 

E  L 

where  £  =  E(AJ)  denotes  the  mean  squared  value  of  the  signal 
amp  1 i tude . 

The  optimum  (minimum  probability  of  error)  receiver 
structure  in  the  presence  of  additive  white  Gaussian  noise 
(however  with  no  jamming  present)  is  the  same  as  that  for 
noncoherent  nonfading  FSK  and  its  structure  is  therefore 
given  by  the  receiver  of  Figure  5. 

B.  ANALYSIS  WITH  NEAR  OPTIMUM  JAMMER 

We  shall  analyze  the  effect  of  Rayleigh  fading  on  an 
incoherent  FSK  receiver  operating  in  the  presence  of  a 
jammer  in  addition  to  the  additive  white  Gaussian  noise. 
Thus,  the  signals  at  the  front  end  of  the  receiver  are 
either 
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Cj  and  Cg.  have  been  defined  in  Equation  3.13.  As  before,  in 
this  equation  SNR  and  JSR  are  defined  by  Equations  2.10  and 
2.11  respectively. 

The  performance  results  to  be  presented  are  a  function 
of  SNR,  JSR,  and  the  number  of  times  the  jammer  sweeps  the 
signal  band.  This  is  discussed  in  greater  detail  in  Chapter 


We  note  that  all  the  terms  in  Equation  3.24  and  Equation 
3.25  take  on  the  form  of  a  SINC  function  times  a  sine  func¬ 
tion.  We  can  readily  show  that  for  even  values  of  m  and  1, 


(nj>C),  -  ( n .  C)  =  0  1  ; 

Using  the  mathematical  forms  given  by  Equations  3.22,  3.23 

and  3.26,  the  performance  of  the  receiver  with  FM  sweep 
jamming  is  obtained  via  the  use  of  Equations  3.8  and  3.9  and 
evaluated  from  the  expression 

(3.27) 

p*  ptH'>  c'  [ 1  -  teW  Jo1  ,JZT& )  oe] 

*  P(H|)  Ci  )  ae] 

*  P(H«>  f1  -  ir JT(/EX.  JcTpT. )  dgj 

*  P(H">  c'  [-g5r  IJ </^7  ,/cTZ' )  d,f] 


where 


cS,r  ?  shp  f  l  ♦  (  p  >  cos9  ,  JSR 

o4.=  2  sw  [  1  .  2/Sp  j„j(  ^  ,  cos(,  .  JSP 

(2,i  -  (3,*  =  2  SNR  JS R  f  p  ) 
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(3.20) 


(n.  ,S), 

J 


O 


T 


£ 


I  v  i3 


r,-„-cc 


Sin  7 c(^K-ftn4i/£) 
7t/2n  K  -t-rn  +  j2/^ 


Sf-n  7t(£7ik-l/,2) 

^(gnk-  -f/fe) 


and 


(n  .,S\ 


ft: 


!*.<  P 


S  i  r)  [aj 1 3  UJ0)  -Ml  u^.  |  T 
-£■(<*/,  t3^-f  no/ 


r  Si n[jr(w -  Mo) T 

2_  —  lA/0)  -j-  T)  uJ* 


which  becomes 


(n.  ,S\ 

j 


Ri,T 


£ 


I  p 


S;^7L(ar\K-im-ite) 
71^2^  K  +  771  -^/£) 


S;n7t(anK-hl/?) 

7^(an  k  +  ^/2) 


, (3.21) 


If  1  is  an  even  integer,  then  Equation  3.20  and  Equation 
3.21  at  most  contain  two  terms  respectively.  These  terms  can 
exist  only  if  the  argument  of  all  SINC  functions  is  zero. 
Therefore  Equation  3.20  and  Equation  3.21  can  be  simplified 
to  yield 


where  n,  =  l/4k,  n^  =  -  (m  ♦  1/2)  /  2k  and 


where 


(  (3  )  -  V  (?>] 

-  (m  -  1/2)  /  2k 


In  view  of  practical  communication  system  constraints, 
the  integer  m  will  typically  be  much  larger  than  1  because  m 
and  1  represent  the  sum  and  difference  of  the  signaling 
frequencies  respectively.  In  particular  we  see  that 


Since  the  FSK  signal  covers  approximately  the  band  (W 0  - 
4  7C/T  ,  W,  +  47C/T)  (assuming  Wt  >  We),  it  is  reasonable  to 
set 

W5  s^-(  W,  +we)  (3.19) 


so  that  the  instantaneous  jammer  frequency  band  (W5  -  (3  W .  , 
ws  +  ^  Wj  )  completely  covers  the  signal  band  provided  that 


ws  - 


=  w. 


47T 

'T 


=  W 


+ 


4 -71 
T 


This  means  that 


pw.  =  d-c  w,  -  w„>  + 

must  be  satisfied.  Since  for  FSK  signaling  we  have  assumed 
that 


(  w,  wG)  T  =  17C  ;  CW,  +  w0)  T  =  m  7L 

where  1  and  m  are  integers,  and  W*T  =  27Ck,  it  is  apparent 
that 


P-  -^[4  *  1/2] 

Note  now  that  the  integer  k  determines  the  number  of  times 
the  jamming  waveform  will  sweep  the  signal  band  in  one  bit 
interval.  Thus  from  Equation  3.18  and  Equation  3.19,  we 
have 


which  becomes 


2T 


r  stn[i(w.-uj,)+r)iv.]T 


■  £  J*(  f*  \  X,  ,  - 

i(U4-W,U  Hu;. 
g“(3  w,  +  7)o/. 

J 
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where  |3  =  and  0  is  a  deterministic  phase  angle. 

The  instantaneous  jammer  frequency  is 


W;(t)  =  V!s  +  H.  cos  W.t 

r  J  j 

so  that  the  instantaneous  jammer  frequency  covers  the  band 


(WS  -  (3W;  , 


w. 


W-  )  ■  Assume  that  W^T  =  2  7C  1  and  W.  T  = 

V  J  sj 


27Ck,  where  1  and  k  are  integers.  In  order  to  determine 
receiver  performance  in  the  presence  of  such  a  jammer,  the 
parameters  SaK  (i,k  =  0,1)  which  are  a  function  of  nj(t), 
must  be  determined.  Before  so  doing,  we  evaluate 


(nj.S>, 


(3.16) 


K  =  0  ,  1 

where  since  the  fixed  phase  $  represents  a  time  delay  only, 
it  is  set  to  zero  for  computational  ease.  Equation  3.16  can 
be  rewritten  as  follows 


(n.,S)K  --M-fl' os  [  Ws  -  "->t  *  P  sin  V. 


■>o\  1  I 

cos  [(Hs+  WK)t  +  p>  sin  W.tJ  } 


dt 


(3.17) 


By  using  the  well-known  Bessel  function  coefficient  expan¬ 
sion  for  each  cosine  term  in  Equation  3.17,  the  following  is 
obtained 


(nj  ,sV 


^T 


y  ja<  ft  ,i 

-CD  r  L  ( -  u/K)  +  mvj 

S?^[(UJs  +U/J  +  71U41T  " 

+  av)  t  n  w. 


(3.18) 
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and  if  the  second  choice  is  used,  we  obtain 


Pe  =  P(H,)  C,  [  1  -  J~(T ('/cu'o  )  de] 

+  p(h,  )  c2  [  jjf  (Jcii:.  )  de] 

+  p(He)  c2  [  1  -  ~:\fi  c/dxo^,;  ,/cJJ,  )  d<j>] 

+  P(H.)  C,  [  -±-  £o  (JC'JJi  JCi.d'x  >  <ty] 


The  effect  of  varying  the  threshold  on  receiver  perform¬ 
ance  for  the  jammer  waveform  of  Equation  3.14  is  discussed 
in  Chapter  5. 


C.  FREQUENCY  MODULATION  SWEEP  JAMMING 


In  certain  situations,  the  need  to  jam  a  certain 
frequency  band  rather  than  a  discrete  set  of  frequencies 
using  tone  jammers  may  arise.  Therefore  an  FM  sweep  jammer 
will  be  proposed,  analyzed,  and  its  effect  on  the  incoherent 
FSK  receiver  investigated  as  a  function  of  the  number  of 
times  the  jammer  sweeps  the  signal  band  during  the  signaling 
interval  [0,T].  The  mathematical  model  used  for  an  FM  sweep 
jammer  .is 


n  ,(t)  = 
J 


in  £  Wst  •+  Kf  J  a  cos  W 


J.t  dt  ]  0  4  t  £  T 


After  carring  out  the  integration  we  obtain 


n .  ( t ) 

j 
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where 


C  =  —  ~  —  1^.  C  =  _ 

'  (I-  oL)\of  ’  Z  y1 

The  effect  of  a  jammer  that  uses  a  single  tone  only  to  jam 
either  'mark'  or  'space'  signals,  can  be  analyzed  by  setting 


(3.14) 


If  the  first  choice  is  used,  we  obtain 


(3.15) 


Pe  =  P(H,)  C,  [ 

1  -  2^1^  >  d8] 

+  P(H,)  C2  | 

2 A  f  9  5  d9] 

+  pch0>  c2  ; 

[  1  "  2tcjC^  ,*3'«  }  d 

l] 

+  P(Hfl)  C, 

[  o  tycr^r  d»j 

whe  re 

f 

<yl\\~  3  Slip  (1  +  2/rH p  cos  6  +  JSR) 
=  0 
2  srn? 

2  SNR  JSR 


Thus  the  probability  of  error  of  the  unmodified  receiver 
(Equation  3.10)  can  be  expressed  in  terms  of  SNR  and  JSR 
only  (defined  by  Equations  2.10  and  2.11)  as  follows 


where 


srnp  f  2  +  2  y  2JSR  cos  0  +JSP) 
o/*e=  SNR  (2-2  /  2 JSR  cos  <j>  +  JSR) 
#,,=  ft,o  =  SNR  JSR 


Receiver  performance  can  now  be  evaluated  as  a  function  of 
SNR  for  fixed  values  of  JSR. 

To  provide  further  insight  into  the  performance  of  the 
modified  incoherent  receiver,  the  effect  of  varying  can 
be  analyzed  via  computation  of  the  probability  of  error  Pe 
from  Equations  3.8  and  3.9.  In  terms  of  SNR  and  JSR,  Pe 
becomes 


=  P(H  , ) 

Pe, 

+  P(H0) 
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=  P(H , ) 
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Observe  that  Equation  3.10  yields  the  performance  of  a 
conventional  incoherent  FSK  receiver  (with  no  channel 
weighting)  in  the  presence  of  jamming.  If  we  now  use  as  the 
jamming  waveform  the  jammer  which  is  optimum  against  a  cohe¬ 
rent  FSK  receiver,  that  is 


n  .(t) 
j 


~  sin  —  (W,  -  VI)  t  cos  (W  ,  + 

=  ft 


W.)  t 


sm  W,t  -  sin 


Wot 


] 


(3.11) 


then  the  terms  SaK  (i,k=0,l)  in  Equations  3.8  and  3.9,  which 
are  a  function  of  the  jammer  waveform  n.(t),  can  be  computed 
as  follows 


S„  =  [(S,,S),  +  (n.,S),]  +  [CS,  ,01,  +(n.,c),_ 
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where  the  dependence  on  6  is  imbedded  in  the  term  S  . 

Following  exactly  the  same  procedure  used  in  obtaining 
the  expression  for  Pe| ,  it  can  be  established  that  the 
expression  for  Pe0,  which  denotes  the  error  probability  when 
H  0  is  assumed ,  that  is,  Prob.(  V  |  >  V  0  |  H  c  )  takes  the  form 


(3.9) 


(1-  c01+ai 


i  Cl  —  dL ) 
vT"  ^  v  i  -  y  )* + (X. 
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where  only  the  term  S  „  is  dependent  on  •  Therefore  the 
total  average  probability  of  error  (Pe).  can  be  obtained  from 
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assuming  that  the  two  hypotheses  are  equally  likely. 


For  the  special  case  =  1/2,  the  performance  of  the 

optimum  receiver  in  the  presence  of  jamming  is  given  by 
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V.  DESCRIPTION  OF  GRAPHICAL  RESULTS 


A.  GENERAL 

In  this  chapter,  the  analytical  results  of  the  previous 
chapters  are  now  presented  via  graphical  means  based  on  the 
derived  mathematical  expressions  for  receiver  probability  of 
error . 

The  plots  presented  display  the  receiver  probability  of 
error  (Pe)  as  a  function  of  SNR  for  the  various  jammer  wave¬ 
forms  previously  considered  for  a  set  of  JSR  values. 

In  each  plot,  the  case  JSR  =  0  has  been  included  in 
order  to  provide  the  basis  for  comparisons  of  the  jammer 
effectiveness  on  the  receiver  performamce  as  it  relates  to 
additive  white  Gaussian  noise  only  interference. 

B.  ASK  (ON  -  OFF  KEYING) 

The  graphical  results  for  the  incoherent  ASK  receiver 
performamce  are  presented  first.  These  plots  correspond  to 
numerical  evaluations  of  Equations  2.5  and  2.7  for  equally 
likely  hypothesis,  that  is,  P(H,)  and  P(He)  are  equal  to 
1/2. 

The  plot  of  Pe  for  ASK  modulation  is  shown  in  Figure  6.7 
as  a  function  of  SNR  for  fixed  values  of  JSR,  using  a  jammer 
as  specified  in  Equation  2.6  .  Figure  6.7  clearly  shows  the 
'break  point’  phenomenon  in  which  if  JSR  increases  beyond  a 
certain  value  (  0.25  in  this  figure),  Pe  increases  with 
increasing  SNR.  From  this  figure,  one  can  observe  that  16.0 
db  of  SNR  is  required  to  obtain  a  Pe  of  ’  -r  without  jammer, 
i.e.,  at  a  JSR  value  of  0.0.  In  comparison,  it  takes  23.5 
db  of  SNR  to  obtain  the  same  Pe  for  a  JSR  value  of  0.1. 
Thus,  in  the  presence  of  a  jammer  with  a  JSR  value  below  the 


break  point,  we  need  a  larger  SNR  in  order  to  obtain  the 
same  performance  level  of  a  receiver  operating  without  a 
jammer  interference.  However,  in  a  jamming  environment,  a 
JSR  value  above  the  break  point  produces  a  Pe  which 
increases  with  increasing  SNR.  In  fact,  there  is  no  value 
of  SNR  that  can  produce  Pe  of  lC_h  for  a  JSR  above  the  break 
point.  For  the  case  of  ASK  modulation,  the  break  point 
occurs  at  a  value  of  JSR  which  is  approximately  0.25,  as 
obtained  from  Equations  2.8  and  2.9. 

Figure  6.8  shows  the  variable  thresholding  effect  on  the 
jamming  situation  with  JSR  =  0.3  beyond  the  break  point  (JSR 
=  0.25).  Instead  of  using  the  fixed  threshold  as  given  by 
Equation  2.5,  the  variation  of  the  receiver  threshold 
obtained  from  Equation  2.5  by  changing  the  value  of  R  ,  can 
reduce  the  jamming  effect  over  a  restricted  range  of  SNR 
values  as  shown.  However,  since  as  shown  in  Equation  2.13 
the  variation  of  the  value  R  does  not  significantly  affect 
the  threshold  value,  the  variation  over  a  wide  range  of 
values  of  R  does  not  result  in  a  significant  change  in  Pe . 

C.  FSK  WITH  TONE  JAMMER 

This  section  presents  graphical  results  pertaining  to 
jamming  effects  on  FSK  modulation  with  a  single  tone  jammer 
acting  against  one  of  the  two  channels  and  a  jammer 
consisting  of  two  different  tones  acting  against  both  chan¬ 
nels  simultaneously. 

Figure  6.9  corresponds  to  the  performance  of  the  optimum 
FSK  receiver  in  which  ^  =  1/2.  Equation  3.10  or  Equation 
3.12  is  used  to  evaluate  performance  with  the  near  optimum 
jammer  specified  in  Equation  3.11  .  This  jammer  waveform 
can  be  thought  of  as  'mark'  and  'space'  channel  jamming. 
Figure  6.9  shows  a  similar  result  to  that  found  in  the  ASK 
case  except  that  the  breakpoint  occurs  at  a  higher  value  of 
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JSR  than  that  found  for  ASK.  This  breakpoint  occurs  at  a 
JSR  somewhere  between  0.5  and  1.0  as  shown  in  this  plot. 
From  this  figure  it  can  be  noted  that  13.5  db  of  SNR  is  need 
to  obtain  a  Pe  of  for  a  JSR  value  of  0.0,  but  the  same 

Pe  is  obtained  by  increasing  the  SNR  to  16.5  db  for  a  JSR  of 
0.1.  This  demonstrates  that  relatively  low  JSR  values 
require  a  significant  SNR  boost  in  order  to  maintain  a 
certain  desired  Pe  value. 

As  shown  in  Figures  6.7  and  6.9,  comparison  of  ASK  and 
FSK  modulation  reveals  that  FSK  is  somewhat  less  vulnerable 
to  jamming.  However  it  must  be  remembered  that  the  jammer 
waveform  n.(t)  used  in  each  case  is  different. 

j 

The  effect  of  the  single  tone  jammer  on  the  optimum  FSK 
receiver  is  presented  in  Figure  6.10  which  corresponds  to 
the  evaluation  of  Equation  3.15  .  The  single  channel 

jamming  on  either  the  'mark'  or  the  'space'  channel  has  the 
same  effect  insofar  as  single  tone  jamming  is  concerned. 
Therefore  the  effect  of  'mark'  channel  jamming  only  is  eval¬ 
uated  and  plotted.  Note  that  in  Figure  6.10,  19.5  db  of  SNR 
is  required  to  obtain  a  Pe  of  ]_.q~s  as  compared  to  an  SNR  of 
24.5  db  in  Figure  6.9  for  a  JSR  value  of  0.3  with  Pe  = 
also . 

As  expected,  Figures  9  and  10  demonstrate  the  fact  that 
single  channel  jamming  is  less  effective  than  simultaneous 
jamming  of  'space'  and  'mark'  channels  with  a  near  optimum 
j  arame r . 

The  effect  of  a  variable  threshold  on  FSK  will  be 
considered  by  changing  the  value  of  q4  in  the  modified 
receiver  shown  in  Figure  6.6  .  For  a  value  of  q/.  other  than 
1  2,  the  simultaneous  jamming  of  'mark'  and  'space'  channel 
results  in  a  compensation  of  the  other  channel  such  that  the 
jamming  effect  remains  the  same  as  in  the  case  of  oC  =  1/2. 
In  other  words,  it  is  difficult  to  reduce  the  near  optimum 
jamming  effect  by  means  of  a  varying  the  threshold. 
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On  the  other  hand,  when  one  channel  jamming  is  applied 
to  the  receiver,  the  jamming  effect  can  be  reduced  by 
adjusting  the  threshold  with  increasing  SNR  as  shown  in 
Figure  6.11,  particularly  for  JSR  of  1.0  and  various  values 
of  .  A  moment's  reflection  will  reveal  that  the  'mark' 
channel  jammer  increases  the  output  power  level  of  the  upper 
envelope  detector  (see  Figure  6.6)  so  that  when  the  'space' 
signal  is  sent,  the  error  increases.  Therefore,  for  this 
type  of  error,  under  the  assumption  that  the  'space'  signal 
has  been  transmitted,  Pe  can  be  reduced  by  lowering  the 
level  of  the  output  of  the  multiplier  by  using  an  appro¬ 
priate  value  of  qC  ■ 

D.  FSK  WITH  FM  JAMMER 

This  section  presents  the  effect  of  an  FM  sweep  jammer 
using  sinusodal  modulation  on  noncoherent  FSK  signaling.  The 
FM  jammer  was  designed  to  sweep  the  bandwidth  occupied  by 
the  signal  several  times  during  a  bit  interval.  Thus  the 
jammer  effectiveness  was  investigated  as  a  function  of  the 
number  of  times  the  jammer  sweeps  over  the  band  of  the 
signal  during  one  bit-time  interval. 

Figure  6.12  shows  the  result  for  one  sweep  of  the  jammer 
per  bit  interval.  Figure  6.13  shows  the  result  of 
increasing  the  sweeping  to  two  sweeps  per  bit  interval. 
These  plots  show  that  the  FM  sweep  jammer  can  be  more  effec¬ 
tive  by  increasing  the  number  of  times  of  jammer  sweeping 
during  a  bit  interval.  This  can  be  expected  from  the  results 
obtained  in  Equation  3.27  .  In  Figure  6.12,  in  order  to 
obtain  a  Pe  of  ,  16.0  db  of  SNR  is  required  for  JSR  value 
of  0.3,  but  in  Figure  6.13  which  corresponds  to  twice  the 
jammer  sweeping,  the  same  Pe  can  be  obtained  by  increasing 
the  SNR  value  to  19.5  db  for  the  same  JSR  value  of  0.3.  In 
comparison  to  the  previous  case  of  FSK  with  tone  jamming, 
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for  the  same  jamming  environment  (i.e.,  JSR  value  of  0.3) 
Figures  6.9  and  6.10  show  that  24.5  db  and  19.5  db  of  SNR 
are  required  respectively  in  order  to  get  the  same  Pe  of 

L  ^ 

These  different  requirements  of  SNR  value  for  various 
jammer  waveforms  show  that  the  FM  sweep  jammer  can  be  effec¬ 
tive  but  in  general  is  not  as  effective  as  the  near  optimum 
jammer.  Note  that  from  a  practical  point  of  view,  the  added 
complexity  of  FM  jammer  waveform  may  make  it  an  unlikely 
candidate  for  a  replacement  of  the  near  optimum  jammer. 
However,  one  advantage  the  FM  sweep  jammer  has  over  the  near 
optimum  jammer  is  that  the  former  can  spread  its  power  over 
a  large  bandwidth  easily  and  therefore  is  more  effective 
than  the  latter  in  the  case  of  lack  of  exact  information  or 
knowledge  about  the  signal  carrier  frequency. 


VI.  CONCLUSIONS 


The  familiar  model  in  which  Gaussian  noise  is  the  total 
interference  is  not  adequate  when  jamming  or  interference 
signals  are  present  in  the  transmission  environment.  This 
thesis  has  analyzed  the  effect  of  various  deterministic 
jammer  waveforms  in  terms  of  probability  of  error  on  binary 
incoherent  receivers  operating  in  the  presence  of  noise. 

From  the  jammer  point  of  view,  the  goal  is  to  cause  the 
maximum  possible  error  to  the  various  receivers  while  making 
efficient  use  of  its  available  power  (i.e.,  with  fixed 
jammer  power).  For  coherent  receivers,  it  was  proved  that 
the  optimum  jammer  waveform  is  made  of  a  deterministic 
signal  proportional  to  the  difference  of  the  binary  signals 
used  to  carry  the  digital  information.  This  thesis  has 
demonstrated  that  those  optimum  jammers  derived  for  coherent 
receivers  perform  their  function  as  near  optimum  jammers 
satisfactorily  against  incoherent  receivers. 

Therefore,  these  near  optimum  jammers  can  be  concluded 
to  be  one  of  the  most  attractive  candidates  for  efficient 
jamming  of  binary  incoherent  communication  systems.  An 
optimum  jammer  has  not  been  derived  or  analyzed  because  the 
complexity  of  the  expression  for  Pe  makes  it  very  difficult 
if  not  impossible  to  derive  the  optimum  jammer  waveform  in 
closed  form.  Other  jamming  waveforms  such  as  single  channel 
jamming  and  FM  sweep  jamming  showed  its  inferiority  in 
comparison  to  near  optimum  jammer  waveforms  and  their  effi¬ 
ciency  can  be  reduced  by  means  of  appropriate  variations  of 
the  receiver  threshold. 
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ASK  WITH  DETERMINISTIC  JAMMING 

PROB.  OF  ERROR  VS  SNR 


Figure  6.7  Near  Optimum  Jammer  for  ASK  Modulation 


ASK  WITH  DETERMINISTIC  JAMMING 

PE  VS  SNR  FOR  JSR  =  0.3 


Figure  6.8  Variable  Threshold  Effect  on  ASK  Modulation 


FSK  WITH  DETERMINISTIC  JAMMING 

PROB.  OF  ERROR  VS  SNR 
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Figure  6.9  Near  Optimum  J.imner  for  FSK  Modulation. 


H  SINGLE  CHANNEL  JAMMING 


le  Channel  Jamming  for  FSK  Modulation 


FSK  WITH  SINGLE  CHANNEL  JAMMING 

_ PE  VS  SNR  FOR  JSR  =  1.0 


Variable  Threshold  Effect  on  Single  Channel  Jamming  for  FSK 


FM  JAMMING 

EEP/CYCLE  (K  =  1) 


Frequency  Modulated  Jammer  for  FSK  Modulat 


FSK  WITH  FM  JAMMING 

TWO  SWEEP/CYCLE  (K  =  2) 
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Figure  6.13  Frequency  Modulated  Jammer  for  FSK  Modulation  (K 


APPENDIX  A 

DIGITAL  COMPUTER  IMPLEMENTATION  OF  THE  MARCUM  Q-  FUNCTION 


The  Marcum  Q- function  occurs  frequently  in  communication 
problems  involving  incoherent  detection  of  signals  with 
single  or  multiple  observations  in  the  presence  of  noise  and 
jamming.  So  it  is  often  necessary  to  compute  values  for  the 
Marcum  Q-function  which  is  defined  by 

,,2.  2  (A.  1) 

Q(  a  ,  8  )  =  V 

h 

-  1  - 


exp ( -  )  IQCa  V)  dV 

f2  2 

V  exp ( -  v  )  I0(a  V)  dV 

2 


where  1^ ( • )  is  the  modified  Bessel  function  of  zero  order. 
It  is  noted  that  the  integrand  f ( V ,  p()  of  Equation  A.l  is 
the  Rician  density  function  which  is  sometimes  called  the 
generalized  Rayleigh  density  function.  The  normalized 
Rician  distribution  for  2  =  1  is  shown  in  Figure  A.l  . 
For  large  values  of  ,  the  Rician  density  function  f(V) 
can  be  approximated  by  the  normalized  Gaussian  density  func¬ 
tion  with  mean  value  of  approximately  ^  .  This  can  be  justi¬ 
fied  as  follows 
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Figure  A.l  Rician  Density  Function. 


by  use  of  the  approximations  IQ(X)  =  ~  ,  and  p(\/  =  V* 

for  large  ^  .  JZK* 

In  order  to  implement  the  computations  of  the  Marcum 


Q-function  using  a  digital  computer,  particular  interest 
must  be  focused  on  values  of  the  function  on  the  tails  of 


the  Rician  density  defined  by  the  integrand  in  Equation  A.l 
Specifically,  small  values  around  the  tails  of  the 


density  must  be  monitored  because  of  the  limitations  of  the 
digital  computer.  By  a  change  of  variables  (X  =  gC.  V)  the 
Marcum  Q-function  defined  above  can  be  expressed  as 


0(  a  ,  S 


(  -  - 

a 


ii2' 


f  (  x )  dx 


(A. 2) 


-x 

where  f(X)  I0(X).  Since  the  digital  computer  can  treat 
the  values  in  the  range  of  y  such  that  -180.218  5  y  5 
174.643  for  the  function  eY  (  for  the  IBM  3033  ),  for 
digital  computer  implementation  of  the  Marcum  Q-function  the 
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integral  limits  need  to  be  adjusted  to  meet  the  conditions 
for  these  acceptable  values  in  the  computer.  Negligible 
areas  around  the  tails  of  the  density  beyond  certain  bounds 
must  be  discarded  without  significantly  increasing  overall 
error.  Therefore  from  Equation  A. 2,  the  integral  limits  can 
be  substituted  with  appropriate  bounds  such  that  exponent 
satisfies  the  condition  given  by 

-  Li >  -J80 

In  other  words,  _ the  variable  X  must  be  located  within  the 

range  -  oi  ^  360  <  X  <  ^  2  *  o(  /  360.  For  a  value  of  X 

beyond  that  range  the  exponential  of  Equation  A. 2  is  too 

- 1  to 

small  (less  than  Q  )  compared  to  other  computed  values 
so  that  the  area  outside  of  the  new  limits  can  be  neglected. 

Let  us  denote  UL  and  LL  the  upper  and  lower  limit 

respectively  defined  by 


UL  = 


a 


LL  =  a ' 


-  a. 


Then  we  can  consider  the  various  situations  case  by  case. 
If  the  lower  limit  LL  is  positive  for  certain  values  of^, 
then  in  the  case  of  (/ft  <  LL,  that  is,  £3  <  -  J 3  6  0 ,  the 
integral  value  can  be  ignored  so  that  the  value  of  Q(^,^  ) 
can  be  assigned  to  be  one.  In  the  case  of  LL  S  £  UL, 
that  is,  pC  -  J 360  S  j3  <  pC  *  J 360,  the  Marcum  Q- function 
can  be  approximated  by 


Q(  ct  ,  S  )  =  1  -  x  exp 

A.L 
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f(x)  dx 


In  the  case  of  ^  >  UL,  that  is,  +  J  360,  the 

Marcum  Q- function  can  be  computed  from 


Q(  a  ,  6  )  =  1  - 
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On  the  other  hand  if  the  lower  limit  LL  is  negative  for  some 
value  of  ,  then  there  are  two  cases  to  be  considered.  In 
the  case  of  £  UL,  that  is,  ^  +  /3  6 0,  the  Marcum 
Q-function  can  be  substituted  with  th_:  value  given  by 


Q(  a  , 

8 

) 

=  1  ■ 

i  r" 

-  — -lx  exp 

r  C  -21  -  a)2  7 

-  _ _ 

2 

Finally  in 

the 

case 

of  *  UL, 

that  is,  ^3 

the  value 

of 

Marcum 

Q-function  can  be  computed 

Q(  a  , 
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=  1 

i  rUL 

-  -  /  x  exp 

,  x  .2  *) 
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-  -  Ci _ 

L  2  J 

f(x)  dx 


f(x)  dx 


Here  it  has  been  assumed  that 

-K 


e  I0(X) 


the  function  f(X)  defined  by 
does  not  impose  the  limitation  of  computation  on 
the  digital  computer  and  that  the  library  functions  for  the 
computation  of  e  Ij(X)  and  its  integration  with  desired 
accuracy  are  available  to  the  user. 
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